N88- 16402 


AUTOMATIC MATHEMATICAL MODELING 
for 

SPACE APPLICATION 
Caroline K. Wang 

Software and Data Management Division 
Information and Electronic Systems Laboratory 
Science and Engineering Directorate 
Marshall Space Flight Center/ NASA 
Huntsville, Alabama 


I. ABSTRACT 

This paper describes a methodology for Automatic Mathematical 
modeling. The major objective is to create a very friendly 
environment for engineers to design, maintain and verify their 
model and also automatically convert the Mathematical model into 
FORTRAN code for conventional computation. 

A demonstration program was designed for modeling the Space 
Shuttle Main Engine Simulation Mathematic Model called 
Propulsion System Automatic Modeling (PSAM) . PSAM is written 
in LISP and MACSYMA and runs on a Symbolics 3670 Lisp Machine. 
PSAM provides a very friendly and well organized environment 
for engineers to build a knowledge base for base equations and 
general information. PSAM contains an initial set of component 
process elements for the Space Shuttle Main Engine Simulation and 
a questionnaire that allows the engineer to answer a set of 
questions to specify a particular model. PSAM is then able to 
automatically generate the model and FORTRAN code. A future 
goal is to download the FORTRAN code to the VAX/VMS system for 
conventional computation . 


II. INTRODUCTION 


Mathematical Modeling for Space Application Simulationpro ject is 
a very complicated process which includes Analysis, Design, and 
the Generation of complex equations. Generally the model will 
require several modification before it will match a real system. 
Historically the modifications have been time consuming and a 
fertile source of error. 
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The use of Artificial Intelligence techniques has shown that this 
process can be simplified. This paper describes a methodology for 
organizing and generating the model automatically. 


III. EXPERT USER INTERFACE 
and 

AUTOMATIC KNOWLEDGE BASE GENERATION 


An Artificial Intelligence system can provide a very friendly work 
environment for engineers. Through PSAM's questionnaire, the 
system collects the necessary information and generates the 
knowledge base automatically. 

The knowledge base includes: (1) An information base in frames and 

(2) Generic equations for the project 


IV. SOFTWARE TOOLS 
and 

AUTOMATIC EQUATION AND CODE GENERATION 


There are a number of tools available for equation derivation . 
and FORTRAN code generation. The one we choose to use is MACSYM 
LISP was the language I used basically for user interface and 
knowledge base generation. I also used LISP for generating 
" MACSYMA" code. The MACSYMA code generates the equations auto- 
matically and at the same time FORTRANcode will also be generated 
automatically and saved into the disk file for future integration. 

The future goal is to integrate the FORTRAN code we generated and 
download to the VAX/VMS system for conventional computation. 


V. MAINTAINING AND VERIFYING 

most difficult problems in software today is the 
and maintenance of existing programs, especially 
proqrams built up our time with many programmers involved. 

This is particularly true for simulation programs. The traditional 
way of modifying simulation programs by rebuilding the model and 
recoding the model had a potential source of many errors. 


One of the 
verification 


234 



The automatic system can eliminates most of these problems. The 
proce-dure become much simpler which makes it easy for the user 
to maintain and modify the model and program directly. 


VI. EXAMPLE 


A demonstration program was designed for modeling the Space Shuttle 
Main Engine Simulation Mathematic Model called Propulsion System 
Automatic Modeling (PSAM) . PSAM is written in LISP and MACSYMA 
and runs on a Symbolics 3670 LISP machine. 

The design goals for PSAM were to develop automatic modeling skills 
for Propulsion System, and other scientific and Engineering applica- 
tions. We used the old Engine Model for an example to study. 

PSAM includes the following features : 

(1) . Friendly user interface. 

(2) . Automatic Knowledge Base generation and 

(3) . Automatic Equation and Coding generation. 

The Space Shuttle Main Engine Simulation model was built up from 
the component process elements and their combination into the sub- 
programs . 

The component process elements are Pump, Hot gas turbine. Hydraulic 
turbine. Turbopump, Combustor, Valve, Incompressible propellant 
flow. Injector volume with priming for start. Hot gas heat transfer 
and Regen cooling flow. The subprograms are Fuel, Oxidizer and Hot 
Gas . 

There are two types of information for a PSAM knowledge base. One 
is the component process elements generic equations and the other 
is the information base for the combination of the Space Shuttle 
Main Engine model subprograms and component process elements . 

The Expert System collects the detailed requirements and generates 


the sets of specific equations for the component process elements 
and subprograms . 

PSAM has the ability to: 

(1) Create or maintain the Knowledge base 

(2) Load different knowledge base 

(3) Automatically generate Equations 

(4) Output generated Equation or FORTRAN code to Disk file or 
option for print out of the Laser printer. 
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interface for the component 


Example for input information user 
process element Pump knowledge base. 



ORIGINAL PAGE IS 
OF POOR QUALXItXi 
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(Deframe Ffpl 

(Class PUMP) 

(Unit "flow-variable") 

(Nomenclature "Low Pressure Fuel Pump Flow variable") 
(Input-parameter Bll DWfdl Sfl)) 

(Deframe Ffp2 

(Class PUMP) 

(Unit " flow-variable" ) 

(Nomenclature "High Pressure Fuel Pump Flow Variable") 
(Input-parameter B18 DWfd2 Sf2)) 

(De frame Fdpl 

(Class PUMP) 

(Unit "flow-variable") 

(Nomenclature "Low Pressure Oxidizer Pump Flow variable") 
(Input-parameter B27 DWmov+DWop3 Sol)) 

(Deframe Fop2 

(Class PUMP) 

(Unit "flow-variable") 

(Nomenclature "High Pressure Oxidizer Pump Flow variable") 
(Input-parameter B39 DWmov+DWotl+DWop3 So2 ) ) 

(De frame Pfdl 

(Class PUMP) 

(Unit "total-pressure") 

(Nomenclature "Low Pressure Fuel Pump Discharge Total pressure") 
(Input-parameter Pfs B12 Sfl Tpfpl)) 

(Deframe Pfd2 

(Class PUMP) 

’ (Unit "total-pressure") 

(Nomenclature "High Pressure Fuel Pump Total Pressure") 
(Input-parameter Pfdl B19 Sf2 Tpfp2)) 

(De frame Podl 

(Class PUMP) 

(Unit "total-pressure") 

(Nomenclature "Low Pressure Oxidizer Pump Total Pressure") 
(Input-parameter Pos B28 S01 Tpopl)) 

(De frame Pod2 

(Class PUMP) 

(Unit "total-pressure") 

(Nomenclature "High Pressure Oxidizer Pump Total Pressure") 
(Input-parameter Podl B40 So2 Tpop2 ) ) 

(De frame Rfpl 

(Class PUMP) 

(Unit "torque") 

(Nomenclature "Low Pressure Fuel Pump "torque"") 

(Input-parameter B13 Sfl Trfpl)) 

(Deframe Rfp2 

(Class PUMP) 

(Unit " torque " ) 

(Nomenclature "High Pressure Fuel Pump "torque"") 
(Input-parameter B20 Sf2 Trfp2)) 

(De frame Ropl 

(Class PUMP) 

(Unit " torque " ) 

(Nomenclature "Low Pressure Oxidizer Pump "torque"") 
(Input-parameter B34 Sol Tropl ) ) 
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PSAM Generic equations : 

•«**** ******PUMP*** ********** "$ 


"Pump Flow Variable"? 

PFV ( F , B , DW , S ) : = 

Block ( F = B * (DW /S) > ; 

"Pump Total Pressure"? 

PTP (PP ,P/B,S,T):= 

Block ( PP = P + B * (S)~2 * T); 

"Pump Torque"? 

PT (R,B,S,T) := 

Block ( R = B * (ST2 * T); 
••♦♦♦♦♦♦♦♦♦♦rpUpBUJE*** ******* "? 


"Turbine Torque"? 


TT(R,B,P,T) := 

Block ( R = B * P * T); 

"Turbine Speed Parameter"? 


TSP(M,B,S,T) := 
Block ( M = B * S / 


(T* ( 1/2 ) ) ) ; 


"Turbine Weight Flowrate"? 


TWF ( DW , B , dP , r ) : = 

Block ( DW = (B * (dP) * r)"(l/2)); 
«********** TURB 0PUMP speed*******"? 


"TurboPump Speed"? 

TPS (S,B, Tq , S0 ) : = 

Block ( S = B * 'Integrate ( Tq,t) 


•«**********VALVE** ******** "? 


"Valve Area"? 

VA(A, Ab,T,Th,Thb) : = 

Block ( A = Ab * (T * Th / Thb))j 


S0) ; 
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"Valve Total Pressure Inlet Total"$ 

VTPIT (P, Pi , RF, DW, Rho) := 

Block ( P = Pi - RF * ( DW ) A 2 /Rho); 

"Valve Total Pressure Inlet Static"$ 

VTPIS(PP,P,DW,A, Rho) : = 

Block ( PP = P - DVT2 / (772 . 8 * A~2 * Rho)); 

"Valve Delta Total Pressure" $ 

VDTP(DP,B,Rhob,Rho,DW,A,Ab,RF) := 

Block ( DP = B * (Rhob/Rho)*(DW/ (A/Ab) )"2 - RF * DW /Rho); 

*********** ********* «»^ 

"Fuel Flow" $ 

FF ( DWf , Bi , P9, P, B j , R, DW0 ) : = 

Block ( DWf = B * 'integrate ( ( (P9 - P - (Bi / R) ) * DWf~2),t) + DW0 ) ; 
"Oxidizer Flow"$ 

OF ( DWo , Bk , Ppos , P f BI , DW0, A, Ab , Bm, DWi ) : = 

Block( DWo = Bk * 'integrate ((Ppos - p - BI * (DWo /(A/Ab))*2 - 
Bm* (DWi) *2 ) ,t) + DW0 ) ; 

"Variable in injector Priming function. Dimensionless "$ 

VIPF (E, Bh , DWo, DWi , E0 ) : = 

Block( E = Bh * 'integrate ((DWo - DWi),t) + E0); 

"Weight Flowrate Injector"$ 

WFI (DWi, DWo, Eo) := 

Block( DWi = DWo * Eo) ; 

****** ****** C ombuster* ********* •*$ 

"Combuster Total Pressure"$ 

CTP (P, Bi , DWI , DW2, B j , DW3, P0 ) : =Block( P = Bi * 'integrate ( ( DWI + DW2 - Bi * 
DW3 ) , t) + p0 ) ; 

"Combuster Fuel Weight Flowrate" $ 


239 



CFWF( F, DW2, DW) : =Block( F = DW2 / (DWl + DW2)); 
"Combuster Temperature"? 

CT(T,R, F,Bk,T9) :=Block ( T = R(F) + Bk * T9); 
«**********COOLING ELEMENT********** "$ 


"Cooling Element Total Pressure"? 

CETP(P,B, DQwl, DQw2,H3, DWi, H,DW,P0):= _ * nw % t \ + p0 \. 

Block( P = B * ' integrated (DQwl + DQw2 + H3 DWl H )» ) ' 

"Cooling Element Specific Enthalpy"? 

CESE(H,B,T) :=Block( H = B * T); 

"Cooling Element Weight Flowrate Main Chamber Heat Exchanger ? 

CEWFM(DW, Bi , P» P5, Bj , R5, DW, DW0) := nw"2^ t) +DW0)- 

Block ( DW = Bi * ' integrate ((P - P5 - (B 3 / R5) * DW 2),t) +DW0J, 

"Cooling Element Density"? 

CED ( R, B , DWi , DW , R0 ) : = . 

Block ( R = B * 'integrate ((DWi - DW),t) + R0); 

"Cooling Element Temperature"? 

CET(T,B,P, R) := 

Block ( T = B * ( P / R))? 

"Cooling Element Heat Transfer Rate Hot Gas Wall"? 

CEHTRHGW ( DQw 1 , B , T , Tw 1 , DWi ) : = Q x 

Block ( DQwl = B * (1.0 + 0.002 * T) * (DWl) 0.8); 

"Cooling Element Heat Transfer Rate Ambient Thrust Chamber ? 

CEHTRAT (DQw2, B, Twl , T, DWi) : = o\. 

Block ( DQw 2 = B * (1.0 + 0.002 * T) * (Tw2 - T) * DWi 0.8), 

"Cooling Element Heat Transfer Rate TC"? 

CEHTR (DQtc, B,Tc, Twl, DWcn) := a 

Block ( DQtc = B * (Tc - Twl) * DWcn 0.8); 

"Cooling Element Hot Gas Wall Temperature ? 

CEHGWT(T, DQtc, DQwl, T0):= * t \ +T 0)* 

Block ( T = B * 'integrate ((DQtc - DQwl) * temp(t),t) + T0J, 
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The Macsyma code that psam generates which will call the Generic 
equation function blocks and also generates the equations and 
FORTRAN code for the certain element. 


fpump( ) : =block(writefile( "als >psam> equation> pump- equation ." ) , 

display(PT (R0P1, B34, SOI, TR0P1 ) ) , 

display(PT(RFP2,B20,SF2,TRFP2) ) , 

display(PT (RFP1 , B13 , SF1,TRFP1 ) ) , 

display(PTP(P0D2, P0D1, B40, S02, TP0P2 ) ) , 

display(PTP(PODl, POS, B28, S01 , TP0P1 ) ) , 

display(PTP(PFD2,PFDl,B19, SF2,TPFP2) ) , 

display(PTP (PFD1, PFS, B 12 , SF1, TPFP1 ) ) , 

display(PFV(F0P2, B39, DWM0V+DW0T1+DW0P3 , S02 ) ) , 

display ( PFV( FDP1 , B27, DWM0V+DW0P3, SOI ) ) , 

display(PFV(FFP2, B18, DWFD2, SF2 ) ) , 

display(PFV(FFPl, Bll , DWFD1 , SF1 ) ) , 

closefile( ) , 

writefile( "al : >psam> fortran>pump . for" ) , 
fortran ( PT ( R0P1 , B34 , S01,TR0P1) ) , 
fortran(PT (RFP2, B20, SF2, TRFP2 ) ) , 
fortran(PT (RFP1 , B13, SF1, TRFP1 ) ) , 
fortran (PTP(POD2, POD1, B40, S02, TPOP2 ) ) , 
fortran ( PTP (POD1, POS, B28, S01 , TPOP1 ) ) , 
fortran (PTP(PFD2, PFD1 , B19, SF2, TPFP2 ) ) , 
fortran ( PTP (PFD1 , PFS, B12, SF1 , TPFP1 ) ) , 
fortran(PFV(FOP2, B39, DWM0V+DW0T1+DW0P3 , S02 ) ) , 
fortran (PFV(FDP1, B27, DWM0V+DW0P3, SOI ) ) , 
fortran ( PFV ( FFP2, B18, DWFD2, SF2 ) ) , 
fortran ( PFV ( FFP1 , B1 1 , DWFD1 , SF 1 ) ) , 
closefile( ) ) ; 


2 

ROP1 = B34 SOI TROP1 
2 

RFP2 = B20 SF2 TRFP2 
2 

RFP1 = B13 SF1 TRFP1 
2 

POD2 = B40 S02 TP0P2 + POD1 


2 

POD1 = B28 S01 TPOP1 + POS 
2 

PFD2 = B19 SF2 TPFP2 + PFD1 


2 

PFD1 = B12 SF1 TPFP1 + PFS 
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FOP 2 = 


FDP1 = 


B39 (DWOT1 + DWOP3 + DWMOV) 
S02 

B27 (DWOP3 + DWMOV) 

SOI 

B18 DWFD2 
SF2 

Bll DWFD1 


FFP2 = 


FFP1 = 


SF1 


ROP1 = B34*S01**2*TR0P1 

RFP2 = B20*SF2**2*TRFP2 

RFP1 = B13*SF1 **2 *TRFP1 

POD2 = B40*SO2**2*TPOP2+POD1 

POD1 = B28*S01**2*TPOP1+POS 

PFD2 = B19*SF2**2*TPFP2+PFD1 

PFD1 = B12*SF1**2*TPFP1+PFS 

FOP2 = B39* (DWOT1+DWOP3+DWMOV) /S02 

FDP1 = B27* (DWOP3+DWMOV ) / SOI 

FFP2 = B18*DWFD2/SF2 

FFP1 = B11*DWFD1/SF1 

The sets or the FORTRAN code will then be integrated into a 
completed compiled SSME program and downloaded to VAX/VMS system 
for conventional computation. 
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